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摘 要 :数值 技术 与 智能 算法 的 飞速 发 展 为 结构 内 部 缺陷 的 识别 提供 了 新 途径 。 本 研究 建立 了 扩 
展 有 限 元 法 (XFEM) 与 误差 反 向 传播 多 层 前 馈 (BP) 神 经 网 络 相 结 合 的 裂纹 反 演 分 析 模 型 。 模 型 通 
过 XFEM 正 向 分 析 获 得 的 测 点 位 移 数 据 训 练 BP 神经 网 络 , 在 此 基础 上 利用 该 网 络 进行 裂纹 反 向 识 
别 。 通 过 2 个 典型 算 例 验 证 了 模型 的 可 行 性 和 精度 ,结果 表明 本 研究 的 方法 能 够 准确 反 演 有 裂纹 的 
几何 信息 ,与 此 同时 还 探讨 了 测 点 布置 方式 及 输入 数据 噪声 等 因素 对 识别 精度 的 影响 。 
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Intelligent crack identification based on XFEM 
and BP neural network 


MAO Xiaomin,ZHANG Huihua,JI Xiaolei ,HAN Shangyu 


(School of Civil Engineering and Architecture ,Nanchang Hangkong University ,330063 Nanchang, China) 


Abstract :The rapid development of numerical technology and intelligent algorithm provides a new way to 
Sidentify the internal defects of structures. In this paper,an inverse analysis model for crack detection is es- 
tablished by combining extended finite element method ( XFEM) and error-back-propagation multilayer 

feedforward (BP) neural network. The BP neural network is trained by the displacement data obtained 
from the forward analysis of XFEM. On this basis , the network is used for the inverse identification of 
cracks. The feasibility and accuracy of the model are verified by two typical examples. The results show 
that the proposed method can accurately retrieve the geometric information of cracks. At the same time ,the 
influence of the layout of measuring points and the input data noise on the identification accuracy is also 
discussed. 
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各 类 工程 结构 因 施工 质量 .外 加 和 荷载 以 及 温度 
变化 等 诸多 因素 的 影响 ,难免 存在 裂纹 和 和 孔洞 等 缺 
陷 。 作 为 代表 性 的 一 类 缺陷 ,裂纹 的 存在 将 导致 结 
构 的 承载 力 下 降 .适用 性 和 耐久 性 退化 ,而 裂纹 的 不 
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断 演 化 甚至 可 能 会 引起 建筑 物 垮塌 等 更 为 严重 的 后 

果 , 危 及 生命 财产 安全 。 因 此 ,对 裂纹 进行 快速 、 准 

确 的 检测 识别 具有 十 分 重要 的 意义 。 
裂纹 识别 属于 典型 的 反 分 析 问 题 。 


反 分 析 是 基 
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于 相关 的 测试 数据 ,通过 建立 有 效 的 反 演 分 析 模 型 ， 
在 不 破坏 结构 本 身 的 前 提 下 识别 出 裂纹 的 位 置 ` 尺 
寸 和 类 型 等 参数 。 目 前 较为 常见 的 裂纹 识别 方法 主 
要 和 包括 图 像 识 别 法 和 无 损 检测 技术 等 ,但 这 些 方法 
存在 着 一 些 局 限 性 , 比如, 图像 识 别 法 易 受 环境 影 
响 ,无损 检测 技术 费用 较为 昂贵 等 。 近 年 来 ,数值 技 
术 与 人 工 智能 的 发 展 为 裂纹 检测 识别 提供 了 一 条 新 
途径 ,此 类 方法 具有 高 效率 ,高 精度 、 低 成 本 和 适用 
性 强 等 特点 。 

基于 数值 技术 和 人 工 智能 算法 进行 裂纹 反 演 包 
括 正 向 分 析 和 目标 极 小 化 两 个 部 分 。 正 向 分 析 主要 
通过 数值 方法 获取 含 裂纹 结构 中 测 点 的 响应 值 (如 
位 移 ) 。 代 表 性 的 数值 方法 包括 有 限 元 法 、 边 界 
元 法 中 .无 网 格 法 中 和 扩展 有 限 元 法 "(extended 人 fi 
nite, element method , XFEM ) 等 。XFEM 模拟 裂纹 问 
题 的 主要 优势 有 :GD 网 格 不 须 与 裂纹 一 致 ,可 在 保持 

网 黎 不 变 的 情况 下 对 不 同 裂纹 构 型 进行 求解 ,因此 
可 高 效率 获取 反 分 析 需 要 的 测 点 数据 ;名 通过 引入 
玛 渤 函数 可 准确 表征 跨越 裂纹 面 时 位 移 的 不 连续 性 
种 绕 尖 的 奇异 性 ,可 在 较 稀 政 的 网 格 上 获得 高 精度 
的 解答 ,因此 能 满足 反 演 类 问题 对 粗 网 格 高 精度 的 
要 驱 。 至 今 ,XFEM 已 被 应 用 于 求解 诸多 裂纹 问 
和 9 目标 极 小 化 主要 通过 智能 优化 算法 来 实现 ， 
常 如 的 智能 优化 算法 有 遗传 算法 '* ( genetic algo- 
rithm, GA) 、 人 工蜂 群 算法 '” (artificial bee colony al- 
g6riihm,ABC) 和 各 类 人 工 神经 网 络 算法 (如 多 层 前 
馈 (BP) 神经 网 络 、 卷 积 神经 网 络 .小 波 神经 网 
络 叶 ' ) 等 。 

QO 近 几 年 来 已 有 许多 学 者 结合 XFEM 与 智能 优化 
算法 进行 缺陷 反 演 。 文 献 [ 12] 提出 了 一 种 基于 
XFEM 和 CA 的 检测 算法 。 文 献 [13] 基 于 XFEM 和 
GA 算法 通过 获取 输入 信号 到 达 测 点 的 时 间 信息 来 
识别 结构 中 的 裂纹 。 此 外 ,文献 [14] 建 立 了 XFEM 
与 GA 相 结合 的 缺陷 反 演 分 析 模 型 。 文 献 [ 15 ] 结合 
XFEM 和 ABC 算法 建立 了 结构 内 部 缺陷 (夹杂 ) 及 
裂纹 的 反 演 分 析 模 型 。 文 献 [16] 基于 XFEM 和 改 
进 ABC 算法 对 结构 中 多 个 缺陷 进行 了 检测 和 量化 。 
文献 [17] 提出 了 一 种 基于 动态 XFEM 和 改进 ABC 
的 多 缺陷 检测 法 。 

如 上 所 述 , 基 于 XFEM 和 智能 优化 算法 (如 GA 
和 ABC) 开展 结构 缺陷 反 演 识别 的 研究 工作 已 有 不 
少 。GA 适用 于 复杂 的 优化 问题 .可 获取 全 局 最 优 解 
且 算 法 独立 于 求解 域 ,但 存在 收敛 速度 慢 .局 部 搜索 
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能 力 差 和 控制 变量 多 等 问题 ;ABC 采用 多 角色 分 工 
和 协同 工作 机 制 ,搜索 精度 较 高 , 鲁 棒 性 较 强 ,但 易 
限于 局 部 最 优 解 且 易 产生 中 间 停 滞 问 题 。 作 为 当前 
应 用 最 为 广泛 的 神经 网 络 模型 之 一 ,BP 神经 网 络 具 
有 强大 的 非 线性 映射 能 力 自学 习 和 自 适应 能 力 . 泛 
化 能 力 .容错 能 力 , 以 及 结构 简单 .工作 状态 稳定 、 易 
于 硬件 实现 等 优点 52 。 考 虑 到 XFEM 模拟 裂纹 问 
题 的 优势 .BP 神经 网 络 的 特点 以 及 尚未 见 两 者 在 裂 
纹 识别 方面 的 协同 研究 ,本 研究 将 XFEM 与 BP 神经 
网 络 算法 进行 有 机 结合 开展 裂纹 反 演 分 析 , 在 验证 模 
型 可 行 性 和 精度 基础 上 ,细致 探讨 测 点 布置 和 输入 数 
据 噪声 等 因素 对 预测 结果 的 影响 ,为 工程 结构 中 裂纹 
的 快速 准确 识别 提供 一 定 的 理论 和 技术 支撑 。 


1 反 分 析 模 型 简介 


对 正 问题 ,一 般 是 在 给 定 材料 类 型 .几何 构 型 和 
边界 条 件 等 情况 下 ,通过 求解 控制 微分 方程 获取 结 
构 的 位 移 、 应 力 和 应 变 等 响应 量 。 而 对 于 反问 题 , 则 
是 通过 已 知 结构 中 某 些 关键 点 ( 即 测 点 ) 的 位 移 等 
响应 量 ,运用 相关 方法 来 反 演 与 该 结构 相关 的 一 些 
未 知 特性 (如 材料 参数 .几何 特性 和 荷载 类 型 等 ) 。 
以 裂纹 识别 问题 为 例 , 其 主要 目的 是 通过 测 点 的 响 
应 量 来 反 推 错 纹 的 构 型 (如 长 度 、 角 度 、 数 量 和 位 置 
等 )。 据 此 ,对 图 1 给 出 的 单一 裂纹 的 识别 问题 , 竺 
反 演 的 参数 向 量 为 

a= {X,Y ,X,Y } (1) 
式 中 :(X 7) 和 (Xe 7) 分 别 为 裂纹 起 点 和 终点 
的 坐标 ,两 裂 尖 均 需 在 结构 体 范围 内 。 


图 1 传单 裂纹 结构 及 测 点 示意 图 


Fig.1 Schematic diagram of structure with a single 


crack and the measuring points 


相应 的 目标 函数 为 
0(a)= 了 du(@a) -wl a 


I 
式 中 :wu 为 由 反 演 算法 得 到 的 测 点 (参见 图 1) 的 响 
应 量 佑 计 值 ;w 为 响应 量 的 真实 值 。 
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据 此 ,裂纹 反 演 分 析 的 主要 目的 是 通过 相关 方 
法 ,在 给 定 的 结构 体 范 围 内 ,获取 最 优 的 参数 向 量 


他 = {证 ，, 产 , ,证 , ,了 ,使 得 式 (2) 最 小 化 。 


2 扩展 有 限 元 法 


2.1 有 裂纹 问题 的 位 移 逼 近 


扩展 有 限 元 法 ”由 美国 西北 大 学 Belytschko 教 
授 在 1999 年 提出 ,是 在 标准 有 限 元 法 (FEM ) 的 基础 
上 发 展 起 来 的 一 种 新 型 数值 方法 。 该 方法 继承 了 
FEM 的 优点 ,但 又 具有 独特 的 优势 ,特别 适用 于 求 
解 如 夹杂 和 裂纹 等 不 连续 问题 。 通 过 在 FEM 的 框 
架 内 引入 表征 裂 面 位 移 不 连续 性 的 跳跃 函数 和 表征 
裂 尖 应 力 场 奇异 性 的 裂 尖 渐 近 基 酉 数 ,对 单元 插值 


丽 数 进行 改进 。XFEM 在 模拟 裂纹 问题 时 不 要 求 网 
格 与 裂纹 面 一 致 ,也 不 需要 使 用 裂 尖 奇异 单元 , 且 在 
较 稀 琉 的 网 格 下 可 获得 高 求解 精度 。 

一 ,对 各 向 同性 材料 中 的 裂纹 问题 ,XFEM 的 位 移 


CN i 
逼近 可 表示 为 
EN mi Na 


> ‘(CY) = > NROu, 十 > NWOH, Ha, 十 
i=1 = 
Na 4 


SNR) DY Bi Hb! (3) 
8=1 gq 


式 击 :X=(X,Y) WN(X),H (XU 和 D(X)(g=1~4) 
分 生 为 与 结 点 iy. 有关 的 单元 形状 函数 、 裂 面 改 进 
函 赂 和 裂 尖 改进 函数 ;N, .Vs 和 Ns 分 别 表 示 单元 内 
的 到 规 结 点 、 裂 面 改 进 结 点 和 裂 尖 改进 结 点 的 数目 ; 
u,b! 分 别 代表 常规 自由 度 、 裂 面 改进 自由 度 和 
裂 尖 改进 自由 度 。 
瓦 (于 ) 和 8( 马 ) 的 具体 形式 为 
1 p(X,0) >0 


Oe) | -1 p(X,0) <0 人 


D; =Vrsin 二 


D? =Vrcos 
2 
8 (5) 
G; =Vrsingsin 7 


Cig =wrsingcos 人 
式 中 :gp (站 ,0) 为 裂 面 水 平 集 函 数 '" ;r 和 9 为 裂 尖 
局 部 极 坐 标 。 

结 点 类 型 的 确定 方式 可 参考 图 2。 具 体 来 说 ， 


Hi 


六 上 VY、 和 人 人 帮 甘 昌 工 
ChinaXiv 器 | F 期 刊 | 第 39 卷 


对 于 被 裂纹 部 分 切割 的 单元 ,这 些 单元 的 结 点 均 为 
裂 尖 改进 结 点 ;对 于 被 裂纹 完全 切割 的 单元 , 除 属 于 
裂 尖 改进 结 点 以 外 的 结 点 均 为 裂 面 改进 结 点 ;除去 
改进 结 点 后 的 剩余 结 点 则 为 常规 结 点 。 


〇 常规 结 点 
图 裂 面 改进 结 点 
@ 裂 尖 改进 结 点 


到 2 XFEM 结 点 类 型 的 确定 
Fig.2 Determination of XFEM node type 


2.2 离散 方程 的 建立 
将 式 (3) 中 的 位 移 近 似 函 数 代 入 虚 功 方程 ,可 
得 到 XFEM 求解 裂纹 问题 的 总 体 离散 方程 , 即 
Kd=F (6) 
式 中 :d 为 结 点 未 知 量 和 矩阵 ;KK 为 整体 刚度 矩阵 ,由 
单元 刚度 矩阵 天 集成 得 到 , 即 
Lp 帮 阿 
天 =| hk ke Ee (7) 
pe pe pe 


式 中 
局 = | (BD)"DBdQ(r,s = ua,b) (8) 
式 中 :D 为 弹性 矩阵 ;B'(r =u,a,b) 的 具体 形式 为 


Nx 0 
B'=, 0 NWN,y 
Ny Nx 
(NH;) x 0 
Br = 0 (NH,) y (9) 
(NHi)y (NH;) x 
(NB!) x 0 
B; = 0 (CVG?) ,|(q=1~4) 
(NDB), (ND') x 


下 为 等 效 结 点 和 荷载 向 量 , 由 各 单元 等 效 结 点 荷 
载 集 合 F" 而 成 ,具体 形式 为 
F=f ff ff (10) 
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式 中 


f° = NAT + [Nfsdn 
11 = NASAT + {NHS,AN CT 


f° = | NGAT + | NOsdn 


式 中 :2, 刀 分 别 表 示 单 元 e 对 应 的 离散 域 .外 力 边 
界 ;f, 和 万 分 别 为 外 力 和 体力 。 


3 ”BP 神经 网 络 


3.1 BP 神经 网 络 的 架构 


__BP 神经 网 络 由 Rumelhart 和 Mccelland 的 科学 
小 组 提出 ,是 一 种 信号 前 向 传播 .误差 反 向 传递 的 多 
层 船 馈 网 络 " 。 它 主要 由 一 个 输入 层 一 个 或 多 个 
隐藏 层 以 及 一 个 输出 层 构成 。 每 一 层 上 的 神经 元 都 
通过 权 值 与 邻 层 上 的 各 神经 元 相互 连接 ,而 同一 层 
的 神经 元 之 间 相 互 独立 。 隐 藏 层 和 输出 层 上 的 各 个 
神 笃 元 都 对 应 于 一 个 激活 函数 和 一 个 阔 值 ;输入 层 
坊 的 神经 元 其 阔 值 为 零 ,其 输出 等 于 输入 。BP 神经 
网 站 中 ,输入 与 输出 的 关系 是 一 个 高 度 非 线性 映射 
基 重 ,如 果 输 入 节点 数 为 mn, 输 出 节点 数 为 9, 则 网 络 
是 炎 n 维 欧 氏 空间 到 g 维 欧 氏 空间 的 映射 。 这 一 映 


S H,, = VI +O +C (12) 
式 中 :及 ,,, .1 和 0, 分别 为 隐藏 层 、 输 入 层 和 输出 
层 钢 神 经 元 个 数 ;C 为 1 ~10 的 常数 。 

图 3 给 出 了 一 个 典型 的 三 层 BP 神经 网 络 ,图 4 
示例 了 对 应 的 隐藏 层 神 经 元 模型 ,此 外 ,输出 层 的 神 
经 元 模型 与 图 4 的 类 似 。 


输入 层 隐藏 层 输出 层 


图 3 BP 神经 网 络 的 结构 模型 
Fig.3 Structure of BP neural network 
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每 个 神经 元 ( 如 神经 元 站 接收 其 他 神经 元 (如 


神经 元 外 的 信息 传递 ,对 神经 元 7 来 说 ,总 输入 为 
b; = A wax; — 0)) (3) 
i=] 


式 中 :wj 为 神经 元 i 和 神经 元 7 的 结合 强度 , 称 为 接 
连 权 值 ;x; 为 输入 信号 ; 0, 为 神经 元 j 的 偏差 ( 阐 
值 ) ;A ) 为 激活 函数 ;b) 为 神经 元 输出 信和 号。 


图 4 ”隐藏 层 的 神经 元 模型 
Fig.4 Neuron model of hidden layer 
由 图 4 和 式 (13) 可 知 ,人 工 神经 元 可 以 看 成 由 
一 组 连接 一 个 加 法 器 .一 个 激活 函数 3 种 基本 元 素 
组 成 。 目 前 用 的 最 多 的 激活 函数 是 单 极 性 Sigmoid 
函数 , 妈 


1 
1+e 


f(x) = (14) 


3.2 ”BP 神经 网 络 的 学 习 方 法 


BP 神 经 网 络 的 训练 过 程 主要 包括 以 下 几 个 
步骤 ” 。 

1) 网 络 初始 化 

中 确定 输入 向 量 :z = [xi ,zo,…,x,] (为 输入 
层 神 经 元 个 数 ) 。 

人 确定 输出 向 量 》 和 期 望 输出 向 量 o, 即 
了 = [7 (9 为 输出 层 神 经 元 个 数 ) ;o = 
[o oo]"。 

确定 隐藏 层 输出 向 量 :b = [2 
为 隐藏 层 神经 元 个 数 ) 。 

初始 化 输入 层 至 隐藏 层 的 连接 权 值 wy = 
[ao yo Was Ws] j=1,2,. ,po 

器 初始 化 隐藏 层 至 输出 层 的 连接 权 值 wh = 
[oa 3 Wars Wi] ,k=1,2,..,go 

2) 隐藏 层 答 出 5 的 计算 

中 计算 隐藏 层 各 神经 元 的 激活 值 y , 即 


3 2 wi — 0;,(] = 1,2,.…,p) (15) 


包 计 算 隐 藏 层 7 单元 的 输出 值 。 将 式 (15 ) 的 激 
活 值 代 入 激活 函数 式 (14) 中 ,可 得 隐藏 层 / 单元 的 
输出 值 为 


,b,] (P 
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,1 = 1,2,.…,p) 
( -Ze)) 
(16) 
式 中 :wj 为 输入 层 至 隐藏 层 的 权 值 ;0 为 隐藏 层 单元 
的 国 值 。 
3) 输 出 层 输出 y; 的 计算 
(计算 输出 层 各 神经 元 的 激活 值 s , 即 


— 0,(k = 1,2,.. “,g) 


计算 输出 层 单元 的 实际 输出 值 六, 即 
1 


p 
(DR) 
1 +e 训 hj 


(17) 


ys = As = ,(k = 1,2,..,9) 
(18) 

式 申 :wi 为 隐藏 层 至 输出 层 的 权 值 ;9, 为 输出 层 单 

元 攀 闵 值 。 

二 4) 判断 算法 是 否 迭 代 结 

@ 判 断 网 络 输出 y 与 期 望 输出 o 之 间 的 误差 是 

从 满足 精度 要 求 。 若 满足 , 则 网 络 迁 代 结束 ; 若 不 满 

感 , 册 利 用 步骤 5) 和 6) 更 新 权 值 和 阔 值 ,然后 通过 

式 N5) ~ (18) 重 新 计算 网 络 输出 值 并 与 期 望 值 进 

往 谍 较 , 依 此 循环 直至 达到 精度 要 求 。 

5 ) 权 值 更 新 

a -BD(1) +BD(1 -1)], 

(ba ,2 ，…… ,有 ) 


bi 


(19) 
WB +1) = (0) + 0 -BD (1) +BD' (1 -1)], 
01,2,…,7) 
(20) 
式 中 :7 为 学 习 速 率 ,有 征 m>0; D(1) = -09J/0w;(?i)， 
D (i) = -9J/9wi(1);B 为 动量 因子 , 且 0<B <1。 
6) 闽 值 更 新 
根据 网 络 输出 六 与 期 望 输出 0 之 间 的 误差 ,分 
别 更 新 0; 和 0;。 


0(t+1) = 0,(1) + nb;(1 一 0b) ,Wt ( OF 一 yx) 
(21) 
0.(t+1)=0(t) + (0 —Y;) (22) 


4 XFEM-BP 神经 网 络 方法 的 实现 


本 研究 结合 XFEM 与 BP 神经 网 络 算法 进行 裂 


应 用 力学 学 
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纹 反 演 分 析 。 该 方法 主要 包括 两 部 分 :第 一 部 分 是 
运用 XFEM 进行 正 向 分 析 以 获取 测 点 的 位 移 啊 应 
量 , 男 一 部 分 则 是 利用 测 点 位 移 数据 对 BP 神经 网 络 
进行 训练 并 据 此 识别 裂纹 几何 参数 。 方 法 的 主要 实 
现 流 程 如 图 5 所 示 。 


其 于 XFEM-BP 神经 网 络 
的 裂纹 智能 识别 


目标 函数 优化 


建立 XFEM 断裂 分 析 方 程 
验证 XFEM 的 精度 


获取 测 点 的 位 移 响应 量 


将 正 分 析 得 到 的 位 移 响 应 
量 及 对 应 的 裂纹 构 型 嘻 入 训练 好 的 BP 神经 : 
BP 神经 网 络 模 型 中 预测 数据 对 网 络 进行 验证 


大 测 得 到 i 


实际 值 进行 比较 计算 两 
对 误差 


修改 相关 参数 


图 5 XFEM-BP 神经 网 络 方法 进行 裂纹 识别 的 流程 图 
Fig.3 Flow chart of crackdetection with XFEM-BP 


neural network method 


S 算 例 分 析 


本 节选 取 2 个 典型 的 裂纹 识别 算 例 ( 单 边 裂 纹 
和 中 心 斜 到 ea i 了 考察 。 为 保证 BP 
神经 网 络 输入 数据 的 可 靠 性 ,首先 对 XFEM 正 向 分 
析 的 精度 进行 验证 ;在 此 基础 上 ,进一步 探讨 测 点 分 
布 情况 和 输入 数据 噪声 对 BP 神经 网 络 预测 结果 的 
影响 。 


5.1 和 矩形 板 内 的 单 边 


如 图 6 所 示 , 和 矩形 板 高 及 =6m, 宽 WW=2m, 板 左 
侧 在 距 板 上 边缘 H/2 处 有 一 长 度 为 a 的 单 边 水 平 裂 
纹 。 板 顶部 受到 ex =1 MPa 的 拉 应 力作 用 , 板 底 边 的 
坚 向 位 移 和 底 边 左 端 的 水 平 位 移 均 为 零 。 杨 氏 模 量 
已 =2.1x102Pa, 泊 松 比 > =0. 3, 按 平面 应 变 进行 
分 析 。 


水 平 裂纹 识别 
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6 单 向 拉 伸 作 用 下 含 单 边 裂纹 的 矩形 板 
Fig.6 A rectangular plate with an edge 
crack under uniaxial tension 

5.1.1 XFEM 精度 验证 

放 对 裂纹 的 正 向 分 析 问 题 , 一 般 通过 校 核 断裂 参 
数 和 的 求解 精度 来 测试 方法 的 可 行 性 。 因 此 ,此 处 主 
学 察 应 力 强度 因子 的 XFEM 仿真 精度 。 根 据 文献 
[283 , 当 HAW=1 且 a/W<0.6 时 , 裂 尖 应 力 强 度 因 


ER 的 计算 公式 为 

A Ki = FovTa (23 ) 
和 P=1.12-0.23| 记 +10.6 [外 ) - 

©O 3 4 

AN 21.71 2 +30. 38 [| (24) 
.二 XxXFEM 模拟 时 ,采用 交互 积分 方法 "站 计算 Ki。 
><1) 收敛 性 测试 


@@ 取 裂纹 长 度 为 0.5m, 分 别 基于 3 种 由 4 结 点 矩 
晶 元 组 成 的 网 格 进行 模拟 (单元 数量 分 别 为 120， 
5 和 780 ,相应 的 离散 域 见 图 7) 。XPEM 仿真 得 到 
的 六, 及 对 应 的 参考 解 列 于 表 1。 可 以 看 出 , 随 着 网 


i aX ;wy 会 作 划 : 
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表 1 不 同 网 格 下 的 Ki (单位 :MPa . mm”) 
Tab.1 Kat different meshes (unit:MPa. m'’”) 


XFEM 解 参考 解 
120 单元 435 单元 780 单元 ( 式 (23)) 
1.772 4 1.8262 1.8563 1.8856 


2) 裂纹 长 度 取 不 同 值 时 的 计算 精度 

裂纹 长 度 分 别 取 0.5、0.6.0.7、…1.2 m, 并 统一 
采用 图 7 (ce) 的 网 格 。 图 8 给 出 了 XFEM 的 结果 及 
相应 的 参考 解 , 可 以 看 出 XFEM 仿真 值 与 参考 解 非 
常 吻合 ,表明 XFEM 模拟 具有 良好 的 精度 。 


9 
8 

全 池 

内 ri 

| 和 二 

加 上 且 5 

po 

友和 

人 目 XFEM 解 
1 一 一 参考 解 
0 
05 0.6 0.7 0.8 0.9 1.0 Ld 2 

裂纹 长 度 /m 


图 8 Ki 随 裂纹 长 度 的 变化 
Fig.8 Variation of K with crack length 

5.1.2 BP 神经 网 络 的 训练 和 预测 结 

1) 网 络 的 训练 

BP 神经 网 络 训练 时 ,采用 测 点 位 移 作 为 输入 数 
据 , 将 裂 尖 4 的 式 坐 标 ( 即 已 ) 作为 输出 数据 。 测 点 
的 位 移 通过 在 裂 尖 位 置 已 知 的 情况 下 采用 XFEM 正 
向 分 析 获 得 。 为 此 ,利用 XFEM 并 基于 图 7(c) 的 网 
格 先后 对 裂纹 长 度 取 0.5 .0.6.0.7、…1.4m 等 10 
种 情况 进行 模拟 。 首 先 按 图 9(a) 的 方式 布置 14 个 
测 点 ,各 测 点 的 坐标 值 列 于 表 2, 因 此 输入 层 和 输出 
层 的 神经 元 数 分 别 为 14 个 和 1 个 ,参照 公式 (12)， 


格 的 加 密 ,XFEM 解 逐 渐 趋 于 参考 解 ,很 好 地 展示 了 
方法 的 收敛 性 。 


隐藏 层 数量 取 为 4。 为 初步 判断 网 络 的 学 习 效 果 , 首 
先 采 用 已 有 训练 数据 对 网 络 进行 测试 ,相关 结果 列 于 
表 3。 可 以 看 出 ,对 10 种 给 定 的 裂纹 , 预测 误差 均 在 
1% 以 内 ,表明 搭建 的 网 络 具有 良好 的 自学 习性 。 

表 2 测 点 坐标 


Tab.2 Coordinates of measuring points 


(a) 120 单 元 (b) 435 单 元 
图 7 XFEM 网 格 
Fig.7 XFEM mesh 
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编号 工 坐 标 了 坐标 编号 。 六 坐标 了 坐标 

1 0.0000 0.7692 8 1.3000 6.0000 

2 0.0000 1.5385 9 2.0000 5.2308 

3 0.0000 2.3077 10 2.0000 4.4615 

4 0.0000 3.6923 11 2.0000 3.6923 

(c) 780 单 元 5 0.0000 4.4615 12 2.0000 2.3077 
6 0.0000 5.2308 13 2.0000 1.5385 

* 0.7000 6.0000 14 2.0000 0.7692 
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于 (四 方案 4: 8 个 测 点 (e) 方案 5: 6 个 测 点 介 方 案 6: 4 个 测 点 
四 到 9 测 点 布置 方案 

Fig.9 Layout of measuring points 
© 8 y 2 P 


CN 表 3 BP 神经 网 络 对 XX 的 训练 输出 结果 
Ss Tab. 3 Training results of BP neural network for X, 


涤 弘 长 度 m 目标 和 训练 输出 X4 ”相对 误差 /% 
(Oo.5 0.5000 0.5004 0.077 
C0.6 0.6000 0.6003 0.042 
-0.7 0.7000 0.7000 0.007 
Ovo.8 0.8000 0.8000 0.006 
0.9 0.9000 0.9001 0.011 
1.0 1.0000 1.0003 0.035 
1.1 1.1000 1.1004 0.032 
1,2 1.2000 1.2000 0.000 
1.3 1.3000 1.2999 0.004 
1.4 1.4000 1.3999 0.008 


2) 网络 的 预测 精度 

型 纹 长 度 先后 取 0.5.0.6.0.7…1.4.1.5m 等 
11 种 情况 ,对 每 种 情形 ,分 别 基于 XFEM 并 按 图 9 
所 示 的 6 种 测 点 布置 方案 获取 对 应 的 位 移 值 。 网 络 
搭建 时 ,采用 4a =0.5 m,0.6m,0.7m,…,1.4m 等 
10 种 情形 对 应 的 结果 作为 训练 数据 ,并 将 据 此 搭建 
好 的 BP 神经 网 络 用 于 裂 尖 坐标 预测 ,即将 4 =1.5 m 
时 的 测 点 位 移 值 作为 网 络 的 输入 数据 以 预测 X。 
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相关 结果 列 于 表 4。 
表 4，BP 神经 网 络 的 预测 输出 和 相对 误差 


Tab.4 Prediction and relative error of BP neural network 


测 点 。 测 点 隐藏 屋 目标 预测 相对 
布置 方案 数量 和 局 XX 误差 /% 
数量 
1 14 4 1.5000 1.4756 1.627 
2 12 4 1.4968 0.213 
3 10 4 1.4897 0.690 
4 8 11 1.5071 0.476 
5 6 3 1.5210 1.400 
6 4 6 1.4661 2.260 


注 :隐藏 层 神经 元 数量 由 式 (12 ) 确定 , 表 中 的 值 对 应 于 预测 精 
度 最 高 的 情形 。 
由 表 4 可 以 看 出 ,在 6 个 测 点 布置 方案 中 ,方案 
2.3 和 4 的 误差 均 小 于 1% ,最 大 误差 出 现在 方案 6， 
但 误差 值 也 不 超过 3% 。 
为 了 测试 BP 网 络 的 稳定 性 ,进一步 对 由 XFEM 
模拟 获得 的 输入 数据 添加 噪声 , 即 
nut =MxreM(T +6) (25 ) 
式 中 :zw 为 BP 神 经 网 络 输入 数据 ( 即 输入 的 测 点 
位 移 值 ) ;uxrey 为 XFEM 正 回 分析 获得 的 测 点 位 移 
值 ;6 为 噪声 的 大 小 (百分比 ) 。 
图 10 给 出 了 6 分 别 取 1% .2% 、5% 和 10% 时 
和 4 预测 的 相对 误差 绝对 值 随 噪 声 值 的 变化 情况 ( 测 
点 按 图 9 方案 2 布置 ,隐藏 层 神经 元 数量 取 4) 。 可 
以 看 出 , 随 着 噪声 值 的 增加 , 裂 尖 位 置 预测 的 相对 误 
差 逐 渐 增 大 。 此 外 ,还 可 发 现 当 噪声 强度 等 级 不 大 
于 5% 时, 网络 预测 的 相对 误差 小 于 3% ,而 当 噪 声 
强度 等 级 为 10% 时 ,相对 误差 小 于 6% ,表明 BP 神 
经 网 络 具 有 很 好 容 噪 性 和 重 棒 性 。 
6 


3 


相对 误差 的 绝对 值 /% 


0 1 2 5 10 
噪声 强度 等 级 /% 
图 10 各 噪声 下 预测 XX 的 相对 误差 绝对 值 
Fig. 10 Absolute value of relative error of 


predicted X, under various noises 
5.2 矩形 板 内 的 中 心 斜 裂纹 识别 


如 图 11 所 示 , 和 矩形 板 高 五 =2.5 m, 宽 WW=2 m 
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板 内 中 心 斜 裂纹 的 长 度 为 a, 与 了 轴 的 夹 角 为 6。 其 
他 条 件 均 与 图 6 相同 。 


中 心 斜 裂纹 的 矩形 板 


A rectangular plate with a central oblique 


图 11 单身 拉 伸 荷载 作用 下 含 

Fig. 11 
crack under uniaxial tension 

XFEM 精度 验证 

一 XFEM 模拟 时 , 取 a =0.5 m, 倾 角 B 依 次 取 为 

22-5°、45° 和 67.5°。 同 样 采用 四 节点 矩形 单元 对 板 

进行 离散 。 图 12 给 出 了 B =45° 的 离散 网 格 ,其 中 包 

念 $500 个 单元 和 546 个 结 点 。 表 5 ~6 分 别 给 出 了 

名 4 应 强度 因 于 和 Ki 的 XFEM 解 及 相应 的 


a 六 ie 


很 高 


表 5 裂 尖 4 的 应 力 强度 因子 太 
Tab.5 Kiof crack tip A 


1/2 
纹 倾 / Ki/(MPa.m’) 
wm . 相对 误差 /% 
《 (OO ) 参考 解 数值 解 
22.5 0.1427 0.1375 3.644 
45 0.4608 0.4660 1.128 
67.5 0.7799 0.7879 1.026 
表 6 裂 尖 4 的 应 力 强度 因子 Ki 
Tab.6 Kof crack tip 4 
1/2 
裂纹 倾 / KI/(MPa.m’) 
相对 误差/% 
(°) 参考 解 数值 解 
22.5 0.3297 0.3251 1.395 
45 0.4475 0.4549 1.654 
67.5 0.3146 0.3161 0.477 


5.2.2 BP 神经 网 络 的 训练 和 预测 结果 
1) 网 络 的 训练 
以 B=45° 为 例 探讨 BP 神经 网 络 的 搭建 问题 。 
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网 络 的 输入 数据 仍 为 XFEM 模拟 得 到 的 测 点 位 移 
值 ,输出 数据 为 裂 尖 4、B 的 坐标 , 即 (X ,YY) 和 (Xs， 
Ys)。XFEM 仿真 时 采用 图 12 的 网 格 , 裂 纹 长 度 分 
别 取 为 0. 5、0. 528 .0. 557、0. 585 .0. 613 、0. 641、 
0. 670 和 0. 698 m,12 个 测 点 的 布置 情况 见 图 12。 
据 此 ,网 络 的 输入 层 .隐藏 层 和 输出 层 的 神经 元 个 数 
分 别 取 为 12、13 和 2。 为 验证 网 络 的 可 靠 性 ,采用 训 


练 数据 进行 裂 尖 坐标 预测 ,相关 结果 见 表 7 和 表 8， 
可 以 看 出 ,对 给 定 的 8 种 裂纹 构 型 ,该 网 络 的 预测 精 
度 都 很 高 。 

(Os) (LS 
(0.0.2.5)， (2.0.2.5) 
(0.0,2.0)4 4 (2.0,2.0) 
(0.0,1.5) 4 4 (2.0,1.5) 
(0.0,1.0)4 4 (2.0,1.0) 
(0.0,0.5) 4 4 (2.0,0.5) 
图 12 XFEM 网 格 和 测 点 布置 (6=45°,a=0.5 my) 
Fig.12 XFEM mesh and measuring points (B =45°,a=0.5 m) 


表 7 BP 神经 网 络 对 裂 尖 4 坐标 的 训练 输出 结果 


Tab.7 Training results of BP neural network for crack tip A 


裂纹 长 度 /m 目标 和 训练 输出 X 目标 ”训练 输出 7 
0.500 0.8232 0.8232 1.0732 1.0732 
0.528 0.8132 0.8133 1.0632 1.0632 
0.557 0.8032 0.8031 1.053 2 1.053 2 
0.585 0.7932 0.7932 1.0432 1.0432 
0.613 0.7832 0.7833 1.0332 1.0332 
0.641 0.7732 0.7732 1.0232 1.0232 
0.670 0.7632 0.7632 1.0132 1.0132 
0.698 0.7532 0.7532 1.0032 1.0032 


2) 网络 的 预测 精度 
基于 已 搭建 好 的 BP 神经 网 络 ,将 a =0.726 m， 
0.755 m,0.783 m 时 由 XFEM 得 到 的 测 点 位 移 作 为 
输入 数据 ,分 别 预测 裂 尖 的 坐标 值 ,相应 的 结果 见 表 
9 ~10。 可 以 看 出 ,BP 神经 网 络 预测 的 最 大 相对 误 
差 为 1. 137% , 表明 搭建 的 网 络 具有 较 好 的 泛 化 
能 力 。 
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表 8 BP 神经 网 络 对 裂 尖 B 坐标 的 训练 输出 结果 


Tab.8 Training results of BP neural network for crack tip B 


裂纹 长 度 /m ”目标 Xs 训练 输出 X。 目标 Ya 训练 输出 Ys 
0.500 1.1768 1.1768 1.4268 1.4268 
0.528 1.1868 1.1867 1.4368 1.4368 
0.557 1.1968 1.1969 1.4468 1.4468 
0. 585 1.2068 1.2068 1.4568 1.4568 
0.613 1.2168 1.2167 1.4668 1.466 8 
0.641 1.2268 1.2268 1.4768 1.4768 
0.670 1.2368 1.2368 1.4868 1.486 8 
0. 698 1.2468 1.2468 1.4968 1.4968 


为 进一步 测试 网 络 的 鲁 棒 性 ,仍然 对 输入 数据 


添加 噪声 。 图 13 ~ 14 给 出 了 噪声 值 依次 取 为 1% 、 
2% 5% 和 10% 时 3 种 不 同 裂纹 长 度 ( 即 a =0.726 m， 
0:955 m,0.783 m) 裂 尖 4 .8 坐标 的 预测 结果 及 相对 
误差 绝对 值 的 变化 情况 。 可 以 看 出 , 随 着 噪声 强度 
等 级 的 增加 网 络 预测 的 误差 也 逐渐 增 大 。 此 外 ,在 
同志 噪声 强度 等 级 下 ,Xi 、3 预测 的 相对 误差 绝对 
依 隔 裂纹 长 度 的 增 大 不 断 增 大 ,而 预测 的 相 
对 误差 绝对 值 随 裂纹 长 度 的 增 大 逐渐 减 小 ;另外 ,还 
醋 殿 发现 裂 尖 4 的 误差 均 大 于 裂 尖 B, 这 可 能 与 约 
瑟 误 件 的 不 对 称 和 裂纹 尺寸 的 影响 等 因素 有 关 。 总 
体 来 说 , 当 噪声 强度 小 于 5% 时 ,预测 结果 的 误差 大 
都 六 10% 以 内 ,表明 搭建 的 网 络 具有 较 好 的 鲁 
性 3 得 与 上 例 相 比 ,此 处 噪声 对 预测 精度 的 影响 更 为 
最 蕃 ,这 很 可 能 与 裂纹 构 型 的 复杂 性 有 关 。 

索 9 BP 神经 网 络 对 裂 尖 4 坐标 的 预测 值 和 相对 误差 


Tab.9 Prediction and relative error of 


BP neural network for crack tip A 


应 用 力 


学 学 报 


一 一 ,预测 
一 e 一 ,预测 
一 一 ,预测 
一 ?一 儿 , 预 泊 
一 9 一 XY, 预测 
一 一 多 预测 


相对 误差 的 绝对 值 /% 


的 相对 误差 绝对 值 ， 
的 相对 误差 绝对 值 ， 
的 相对 误差 绝对 值 ， 
的 相对 误差 绝对 值 ， 
的 相对 误差 绝对 值 ， 
的 相对 误差 绝对 值 ， 


ChinaXiVv 合 ; 


a=0.726 m 
a=0.755 m 
a=0.783 m 
a=0.726 m 
a=0.755 m 
a=0.783 m 
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噪声 强度 等 级 /% 
不 同 噪声 下 X Xs 预测 的 相对 误差 绝对 值 
Absolute value of relative error of X and X, 


prediction under different noises 


一 -了 ,预测 的 相对 误差 绝对 值 ， 
一 e 一 7 预测 的 相对 误差 绝对 值 ， 
预测 的 相对 误差 绝对 值 ， 
-了 ,预测 的 相对 误差 绝对 值 ， 
一 了 预测 的 相对 误差 绝对 值 ， 
一 一 了 ,预测 的 相对 误差 绝对 值 ， 


a=0.726 m 
a=0.755 m 
a=0.783 m 
a=0.726 m 
a=0.755 m 
a=0.783 m 


裂纹 长 度 / 目标 预测 相对 目标 预测 相对 
m Xa X14 ”误差 /% Ya J ”误差 /% 
0.726 0.7432 0.7372 0.806 0.9932 0.9868 0.648 
0.755 0.7332 0.7312 0.268 0.9832 0.9798 0.341 
0.783 0.7232 0.7298 0.910 0.9732 0.9843 1.137 


BP 神经 网 络 对 有 裂 尖 坐标 有 的 预测 值 和 相对 误差 
Tab. 10 Prediction and relative error of 


BP neural network for crack tip B 


裂纹 长 度 / 目标 预测 相对 目标 预测 相对 
Xp Xp ”误差 /% Ys Ys ”误差 /% 
0.726 1.2568 1.2628 0.476 1.5068 1.5132 0.425 
0.755 1.2668 1.2688 0.154 1.5168 1.5201 0.216 
0.783 1.2768 1.2702 0.514 1.5268 1.5159 0.714 
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绝对 值 /% 


相对 误差 的 
S 


0 1 2 5 10 
噪声 强度 等 级 /% 
到 14 不 同 噪声 下 站 、Y 预 测 的 相对 误差 绝对 值 


Fig.14 Absolute value of relative error of Y, and Y, 


prediction under different noises 


6 结 论 


通过 引入 裂 面 和 裂 尖 改进 函数 ,扩展 有 限 元 法 
可 使 用 较 稀 玻 的 固定 网 格 对 不 断 变化 的 裂纹 构 型 进 
行 高 精度 模拟 。BP 神经 网 络 优化 算法 则 具有 较 强 
的 非 线 性 映射 能 力 、 自 学 习 和 自 适 应 能 力 \ 泛 化 能 
以 及 容错 能 力 等 优点 。 本 研究 结合 两 者 的 优点 , 建 
立 了 用 于 结构 内 部 裂纹 反 演 分 析 的 XFEM-BP 神经 
网 络 方法 。 在 给 出 XFEM 和 BP 神经 网 络 基本 原理 
的 基础 上 ,重点 阐述 了 两 者 的 结合 策略 和 实现 流程 ， 
最 后 通过 两 个 典型 算 例 对 其 精度 和 和 鲁 棒 性 等 进行 了 
验证 ,并 探讨 了 网 格 密 度 \ 测 点 布置 和 输入 数据 噪声 
等 对 精度 的 影响 。 结 果 表 明 ,XFEM 正 向 分 析 和 BP 
神经 网 络 的 裂纹 识别 精度 均 较 好 ,此 外 搭建 的 BP 神 
经 网 络 还 具有 较 好 的 容 噪 性 能 。 

需要 说 明 的 是 ,虽然 具有 诸多 优点 ,BP 神经 网 
络 也 存在 收敛 速度 慢 和 易 陷 和 人 局 部 极 小 值 等 缺点 ， 
为 提升 BP 的 性 能 ,可 进一步 引入 其 他 算法 (如 GA) 
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对 其 进行 改进 ,相关 工作 正在 进行 之 中 。 
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